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ABSTRACT 

We estimate the power spectrum of mass density fluctuations from peculiar 
velocities of galaxies by applying an improved maximum-likelihood technique to 
the new all-sky SFI catalog. Parametric models are used for the power spectrum 
and the errors, and the free parameters are determined by assuming Gaussian 
velocity fields and errors and maximizing the probability of the data given the 
model. It has been applied to generalized CDM models with and without COBE 
normalization. The method has been carefully tested using artificial SFI catalogs. 
The most likely distance errors are found to be similar to the original error 
estimates in the SFI data. The general result that is not very sensitive to the prior 
model used is a relatively high amplitude of the power spectrum. For example, at 
k = 0.1 /iMpc-i we find P{k)n^-^ = (4.4±1.7) x 10^ {h'^Mpc)^. An integral over 
the power spectrum yields ag^^'^ = 0.82 ± 0.12. Model-dependent constraints 
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on the cosmological parameters are obtained for families of CDM models. For 
example, for COBE-normalized ACDM models (scalar fluctuations only), the 
maximum-likelihood result can be approximated by f2 hgo^'^ = 0.58±0.11. The 
formal random errors quoted correspond to the 90% confidence level. The total 
uncertainty, including systematic errors associated with nonlinear effects, may 
be larger by a factor of ~ 2. These results are in agreement with an application 
of a similar method to other data (Mark III). 

Subject headings: Cosmology: observations — cosmology: theory — dark matter 
— galaxies: clustering — galaxies: distances and redshifts — large-scale structure 
of universe 
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1. INTRODUCTION 

In the standard picture of cosmology, structure originated from small-amplitude density 
fluctuations that were amplified by gravitational instability. These initial fluctuations are 
assumed to have a Gaussian probability distribution, fully characterized by their power 
spectrum (PS). On large scales, the fluctuations are expected to be linear even at late 
times, still characterized by the initial PS. Thus, the PS is a very useful statistics for large 
scale-structure. We focus on the PS rather than the correlation function {e.g., Gorski et al. 
1989) because the PS distinguishes more clearly between the processes that affect structure 
formation on different scales. It also has the advantage of being less sensitive to assumptions 
regarding the mean density. 

The PS has been estimated from several redshift surveys of galaxies (see reviews by 
Strauss & Willick 1995; Strauss 1998). However, the distribution of galaxies does not 
necessarily provide a direct measurement of the underlying mass distribution; the PS 
estimated from redshift surveys is contaminated by unknown "galaxy biasing" . Additional 
complications arise from redshift distortions, triple-value zones and nonlinearities of the 
density field. This also complicates a direct comparison of the correlation function derived 
from galaxy density fluctuations with similar quantities derived from peculiar velocity 
measurements. Therefore, it is advantageous to estimate the mass PS directly from purely 
dynamical data. Another advantage of velocity over density data is that they probe the 
density field on scales larger than the sample itself, and that they are subject to weaker 
nonlinear effects. It is therefore easier to obtain an approximation for the initial PS from 
the current velocity PS than from the current density PS. 

Direct estimation of the PS from reconstructed velocity or density fields is complicated 
by the need to correct for the effects of large noise, smoothing, and finite and nonuniform 
sampling {e.g., Kolatt & Dekel 1997). On the other hand, the likelihood analysis of peculiar 
velocities, such as the one applied here, provides an appealing method for estimating the 
mass PS since it is a straightforward statistic acting on the 'raw' data, without the need for 
processing such as binning, smoothing, or applying a full POTENT reconstruction (Dekel, 
Bertschinger & Faber 1990). It takes into account the measurement errors and finite discrete 
sampling, and it utilizes much of the information content of the data. The simplifying 
assumptions made in our main analysis are that the peculiar velocities are drawn from a 
Gaussian random field, that the velocity correlations can be derived from the density PS using 
linear theory, and that the errors in the measurements are Gaussian. Other limitations of the 
method are the need to assume some parametric functional form for the PS, with a possible 
sensitivity of the results to the choice of this model, and the fact that the likelihood analysis 
provides only relative likelihood of the different models, not an absolute goodness-of-fit. 
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PS estimates using likelihood analysis (Zaroubi et al. 1997) has been obtained from 
the "Mark III" catalog of peculiar velocities (Willick et al. 1997a), yielding relatively high 
values for the PS, in agreement with the direct estimates from the "POTENT" reconstruction 
(Kolatt & Dekel 1997). This result is still associated with large uncertainties because the 
sampling of the data is sparse and nonuniform, because the merging of data from several 
sources is nontrivial, and because the distance errors in peculiar-velocity data are relatively 
large. Furthermore, the uncertainty in the assumed distance errors always propagates into 
an uncertainty in the resultant PS because the errors add in quadrature to the PS. It is 
therefore important to analyze new data of certain improved qualities and to pay special 
attention to the error estimates. 

The data analyzed in the present paper are based on the new SFI catalog of peculiar 
velocities of galaxies (Haynes et al. 1999, Wegner et al. 1999), containing about 1300 field 
spiral galaxies with TuUy-Fisher (TF) distances. Most of the measurements in the SFI 
catalog arc new. Data taken from the literature which are included in the catalog, mostly 
those by Mathewson, Ford & Buchhorn (1992), have been recalibrated to match the new 
observations both for magnitude and line width scale. This procedure should minimize the 
effects of combining different datasets, effects of significant concern in Mark 111 (e.g., Willick 
& Strauss 1998). The SFI catalog, though sparser than Mark III in certain places, covers 
more uniformly the volume out to 70h~^Mpc. 

The distances in the SFI catalog have been estimated using a linear TF relation derived 
from a matching cluster sample (Giovanelli et al. 1997a, 1997b; SCI). Possible deviations 
from the standard, linear TF relation were ignored since no clear evidence for such deviations 
was detected in the data. In addition, the sensitivity to such an effect is small because of 
the the selection criteria of the SFI catalog. 

The crucial issue of error estimate is addressed in two ways. First, the fact that the SFI 
field sample is matched by the SCI cluster sample of similar size allows a careful investigation 
of the observational and internal scatter of the TF distances which provides a good a priori 
estimate of the errors. These errors are adjusted for an assumed difference in the scatter 
between field and cluster galaxies. An additional adjustment of the scatter is due to our 
bias-correction procedure. 

A second and independent approach to estimate the errors is to include them as an 
extra parameter in the likelihood analysis so that it also determines the maximum-likelihood 
values for the errors. In that approach, we use a parametric model for the errors which 
builds upon the original estimates of width-dependent errors. 

We address here the mass-density power spectrum as derived from peculiar- velocity data. 
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with or without Cosmic Microwave Background (CMB) fluctuation data, but independent 
of the distribution of galaxies in redshift space. We thus determine the quantity P{k)Q^ '^ 
(where P{k) is the density power spectrum and Q is the cosmological density parameter), 
while we are free of assumptions regarding the "biasing" relation between galaxies and mass. 
We can therefore measure a purely dynamical parameter such as ag = cxgfi^'^ (where ag is 
the rms mass-density fluctuation in top-hat spheres of radius 8h~^Mpc). When assuming a 
priori a parametric functional form for the mass PS, e.g., based on a generic CDM model, 
we can in fact determine a combination of dynamical parameters such as Q and the power 
index n. 

Investigations involving galaxy redshift surveys commonly measure a different parameter 
that does involve galaxy biasing, (3 = Q^'^/b (where b is the biasing parameter). The 
parameters ag and (3 (at Sh^^Mpc) are related via agg, referring to the rms fluctuation 
in the galaxy number density. A number of measurements of (3 have been carried out so far, 
either based on redshift distortions of the IRAS 1.2 Jy redshift survey (Fisher et al. 1995) 
or based on comparisons of this redshift survey and the peculiar-velocity data. Most recent 
velocity-velocity comparisons found values for (3 in the range of 0.5 — 0.7 (Davis, Nusser & 
Willick 1996; Willick et al. 1997b; da Costa et al. 1998; Kashlinsky 1998; Willick & Strauss 
1998), while density-density comparisons have lead to values as high as 0.9 (e.g. Sigad et al. 
1998). A determination of dg from the SFI data may help to clarify the situation. 

In § 1^ we describe the data and our method for correcting Malmquist bias. In § |^ we 
present the method of analysis and the parametric models used as priors. The method is 
tested using mock catalogs in § ^ The estimated power spectra and the constraints on the 
cosmological parameters are presented in § |^. The robustness of the results is addressed in 
§ ^. We discuss our results and conclude in § ^ 

2. DATA 

2.1. Sample and Distance Errors 

The SFI sample is based on a wide-angle survey of Sbc-Sc galaxies with I-band TF 
distances, covering declinations 5 > —45° and galactic latitudes b > 10°. The galaxy selection 
criteria depend on redshift in order to ensure dense sampling at large distances; the catalog 
consists of three zones of different diameter limits and redshift limits. This data set was 
complemented south of 6 = —45° with galaxies drawn from the Mathewson et al. (1992) 
survey, carefully converted to the same system of magnitude and line-width, and with the 
same set of corrections and selection criteria applied to the whole sample (Giovanelli et al. 
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1997a, 1997b). The combined sample comprises of about 1300 field galaxies, extending out to 
7500 kms~^ in redshift, and quite isotropically covering the whole sky except of the Galactic 
zone of avoidance. 

Accurate estimation of the uncertainty A in the distance are important both for the 
bias correction (see § |2]^) and for determination of the PS (see § The uncertainties are 



derived from the estimate of the scatter in the observed TF relation. We take advantage of 
the fact that the SFI sample is matched by a similar cluster sample (SCI, Giovanelli et al. 
1997a, 1997b). The line- width dependent scatter of this cluster sample is well determined. 
Since SCI was observed using the same observational procedures as most galaxies in SFI, the 
distance estimates in both samples should suffer from similar observational uncertainties. 
However, it is less clear whether the intrinsic scatter of the TF relation is the same for the 
field and cluster samples. We have parameterized the total scatter in the SFI sample by 
using the SCI observed scatter and adding an additional intrinsic scatter for field galaxies in 
quadrature. Such a higher scatter for field galaxies has consistently been found by a number 
of authors (e.g. Bothun & Mould 1987; Freudling, Martel & Haynes 1991). We estimated 
the total scatter of SFI by taking advantage of the distance dependence of biases in the 
inferred distances. These biases for field galaxies are large at high distances and dominate 
the raw measured peculiar velocities. The exact behavior at large distance depends on the 
assumed amplitude of the scatter (see Freudling et al. 1995). With the aid of mock samples, 
the observed distance dependence of the average measured peculiar velocity was used to infer 
the intrinsic scatter for the SFI sample. 

The resulting errors are estimated to be in the range 15 — 20%, and increasing with 
decreasing line-width w. Following da Costa et al. (1996), a small fraction (~ 7%) of galaxies 
with small line-width (logw < 2.25) has been discarded because of the unreliability of the 
TF relation and its scatter at such line-widths. A detailed account on the sample selection, 
error estimates, and the procedure of combining the two datasets can be found in Wegner 
et al. (1999) and Haynes et al. (1999). The SCI sample of ~ 500 galaxies within 24 clusters 
was used for calibrating the TF relation and in estimating the scatter properties, but the 
peculiar velocities of these clusters themselves are not used in this work as they require a 
different treatment. 



2.2. Bias Correction 

It is crucial to properly correct the data for systematic biases, such as those arising 
from the coupling between the random distance errors, the geometry of space and the 
inhomogeneities in the underlying distribution of galaxies, and certain aspects of the sample 
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selection. Due to the complexity of the selection criteria and the TF distance errors in the SFI 
data, the bias correction could not be properly estimated using the standard simple analytic 
expression. In earlier papers of the SFI series, the bias was estimated using a numerical 
Monte-Carlo approach in which the selection criteria were mimicked in detail (Freudling 
et al. 1995). Here, we replace it with a simpler semi-analytic estimate of the bias, which 
incorporates the relevant selection criteria. 

The bias-correction method will be described in detail by Eldar et al. (1998). Here we 
mention only the basic features of the method. Given a galaxy with a TF inferred distance 
d and a line-width rj = \ogw — 2.5, the Malmquist-corrected distance is adopted to be the 
conditional expectation value of the true distance, r, 

!,-drPir,d,r,) ' 

where P{r,d,ri) is the joint probability distribution in the catalog {e.g., Strauss & Willick 
1995). The line- width is explicitly included to ensure that the correction holds when the 
selection criteria depend on i]. 

This joint distribution is derived from several input quantities. One is the underlying 
spatial number density of galaxies, n(r), which is taken from a self-consistent real-space 
reconstruction from the IRAS 1.2 Jy redshift survey (as described in Sigad et al. 1998). 
Another input is the distribution of galaxy diameters, $(-D). One also needs as input the 
conditional probability P{ri, B, I\D), that a galaxy with a given D will have a line-width rj 
and absolute magnitudes B and /. We adopted the same distribution functions as those used 
in Freudling et al. (1995). Taking into account the selection in angular diameter, a = D/r, 
and the apparent blue magnitude limit mB,max, one obtains 

P{r,d,r]) (X r'^n{r) J da S{a\r) ^{ar) r exp f — j P{r,d,T],ar\mB,mii^). (2) 

The selection function of angular diameters at a given true distance, S{a\r), is derived from 
the corresponding selection function in redshift space, Sz{a\z), via 

/oo 
dz S,{a\z) P{z\r), (3) 
-oo 

where P{z\r) is based on the model peculiar-velocity field. The joint distribution 
P{r,d,7],ar\mB,ma.x) is based on a combination of the diameter-magnitude relation, the 
correlation between B and I Magnitudes, the rj — B relation, and the B-magnitude limit 
in the selection of galaxies. 

This bias correction scheme was tested, and its details were refined, using carefully 
constructed mock catalogs (presented below, § [4.1| ). We also tried several variants of the 
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procedure to correct the real data for biases. The results of the power-spectrum analysis 
turn out to be fairly insensitive to the specifics of the bias correction scheme. In particular, 
for the underlying galaxy-density field that enters the correction via n(r), we tried replacing 
the IRAS field (Sigad et al. 1998) with a linear reconstruction of a combination of IRAS 
and optical data (Freudling, da Costa & Pellegrini 1994), and found negligible effects on the 
results of the likelihood analysis. 

The estimated errors in the observed TF relation can be directly translated into 
a distance uncertainty for each galaxy prior to the correction for biases. However, 
the correction for biases changes the properties of the scatter for a given location in 
estimated-distance space, which leads to a different uncertainty in the bias-corrected distance 
estimate. The semi-analytic approach is used also for a re-evaluation of the distance errors 
after the bias correction. We find that the bias correction acts towards slightly decreasing 
(by ~ 10%) the average error, because of the additional information incorporated by the 
selection effects and the underlying density field used in the bias correction. The validity of 
this approach is verified using the mock catalogs (which also shows that the distribution of 
distance errors after the bias correction closely resembles a Gaussian distribution, Eldar et 
al. 1998). An independent verification of the magnitude of errors within the framework of 



the likelihood analysis is described in § In what follows, we refer to these errors as our 
'original' error estimates. 



3. METHOD 

3.1. Likelihood Analysis 

The goal of this paper is to estimate the power spectrum of mass density fluctuations 
from peculiar velocities, by finding maximum likelihood values for parameters of assumed 
model power-spectra. Again, the underlying assumptions are that the velocities and their 
errors are Gaussian, and that the velocity correlations can be derived from the density 
PS using linear theory. The assumption regarding the Gaussianity of the velocity field is 
supported by simulations which show that it is Gaussian well into the quasi-linear regime 
(Kofman et al. 1994). This is farther verified for our data set by the fact that the distribution 
of observed hi{z/d) closely resembles a Normal distribution. The validity of the second 



assumption is discussed later in § H.3.1 . The likelihood analysis method is described in 
Zaroubi et al. (1997; see also Kaiser 1988; Jaffe & Kaiser 1994). Here we summarize the 
main ideas, the underlying assumptions, and the specific application to peculiar velocities. 
Given a data set d, our objective is to estimate the most likely model m. Using Bayes 
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theorem 

and assuming a uniform prior V{m), this can be turned to maximizing the hkehhood 
function, the probabiUty of obtaining the data given the model, £ = V{d\m), as a function 
of the assumed model parameters. 

Under the assumption that both the underlying velocities and the observational errors 
are independent Gaussian random fields, the likelihood function can be written in the 
following form 

C = [(27r)^det(i?)]-V2exp (^-lY^u^R^u^j . (5) 

This is simply the corresponding multivariate Gaussian distribution, where {ui}^i is the 
set of N observed pecuhar velocities at locations {rj}, and R is their correlation matrix. 
Expressing each data point as the sum of the actual signal and the observational error 
Ui — Si + Ci, the elements in the correlation matrix have two contributions 

Rij =< UiUj >=< SiSj > + < CiCj >= Sij + Ci^Sij. (6) 

The first term is the correlation of the signal, that is calculated from theory. The second term 
is the contribution of the distance errors, which are assumed to be uncorrelated. This should 
be true for the observational errors and the intrinsic scatter of the TF relation. We tested 
the impact of uncertainties in the bias correction, which might lead to correlated errors, by 
varying parameters of our bias model within the expected uncertainties. The changes in the 
results reported below are negligible compared to other systematic and random errors. For a 
given PS, the signal terms are calculated using their relation to the parallel and perpendicular 
velocity correlation functions, and ^^j., 

Sij = ^I'±(r) sin 9i sin 9j + ^ty (r) cos 9i cos 9j , (7) 

where r = |r| = |rj — r^l and the angles are defined by = • f (Gorski 1988; Groth, 
Juszkiewicz & Ostriker 1989). In linear theory, each of these can be calculated from the PS, 

*±,||W = ^|^| P{k)K^,\\{kr)dk, (8) 

where K±{x) — ji{x)/x and K\\{x) — jo — 2ji{x)/x, with ji{x) the spherical Bessel function 
of order /. The cosmological fl dependence enters as usual in linear theory via /(Jl) — 
and Hq is the Hubble constant. 



The likelihood analysis is performed by assuming some parametric functional form for 
the PS. For each assumed PS, the correlation matrix R is obtained and used to calculate 
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the likelihood function (eq. 0). Exploring the chosen parameter space, we find the PS 
parameters for which the likelihood is maximized. (Note that since the model parameters 
appear also in the normalizing factor of the likelihood function, through i?, maximizing the 
likelihood is not equivalent to minimizing the x^O The main computational effort is the 
calculation and inversion of the correlation matrix R in each evaluation of the likelihood. It 
is an X matrix, where the number of data points N is typically more than 1000. 

Since the input data are pecuhar velocities, the method essentially measures the 
combination /(f2)^P(/c), and not directly the mass-density P{k) by itself. This degeneracy 
between Vt and the PS can be broken when f2 enters explicitly into the functional form 
characterizing the PS shape, as in CDM models (§ |3.2| ). 

Confidence levels are estimated by approximating — 21n£ as a distribution with 
respect to the model parameters. The likelihood analysis provides only relative likelihoods 
of different models. An absolute measure of goodness-of-fit can be provided, for example, 
by the value of the obtained with the parameter values associated with the maximum 
likelihood. A per degree of freedom of about unity would indicate that the model provides 
a good statistical description of the data. 



3.2. Power Spectrum Models 



In order to perform the likelihood analysis, a specific parametric form for the PS is 
needed. For the main analysis of the paper, we use families of generalized CDM models 
normalized by the COBE 4-year data. The general form of these models is 



P{k) = AcoBE{n, n, A) T\n, Qb, h; k) ¥' 



(9) 



where A is the normalization factor and T{k) is the CDM transfer function proposed by 
Sugiyama (1995, a slight modification of Bardeen et al. 1986): 



T{k) 



ln(l + 2.34g)) 



2.34g 



1 + 3.89g + (16.1g)2 + (5.46g)=^ + (6.71g)^ 



-1/4 



nh exp(-f]b - {2hf^^nb/n){hMpc'^] 



(10) 
(11) 



These models include open universes with no cosmological constant, fiat models with a 
cosmological constant (fi-t-^A = 1), and tilted models with a large-scale power index n that 
can be different from unity. The latter may include tensor fluctuations with tensor to scalar 
ratio of quadrupole moments of T/S = 7(1 — n). The free parameters in the CDM models 
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are fl, h and n. In all cases the baryonic density is set to be Qb = 0.024/i ^ {e.g., Tytler, 
Fan & Buries 1996). For each model, the amplitude A is fixed by the COBE 4-year data. 

We followed the COBE normalization adopted in Zaroubi et al. (1997), who used the 
COBE DMR data (Hinshaw et al. 1996) to set the PS amplitude calculated by different 
authors (Gorski et al. 1995; Sugiyama 1995; White & Bunn 1995) for various cosmological 
CDM-like models. The calculation of Sugiyama (1995) was used as a reference. For models 
not studied by him the other results were used, after matching them to Sugiyama's in the 
commonly studied models. For a summary of the COBE normalization results see Gorski et 
al. (1998). 

In addition, we use a different parameterization of the same power spectra, namely 

P{k)=AkT\k), T{k) = (l + [aklT + {hk/Tf/^ + {ck/TfYy^'\ (12) 

with a = 6.4h-^Mpc, b = S.Oh'^Mpc, c = l.Th'^Mpc and u = 1.13 {e.g., Efstathiou, 
Bond & White 1992). In the context of the CDM model, F has a specific cosmological 
interpretation, F = Qh. Below, however, we use equation ([T^ ) as a generic form with 
limiting logarithmic slopes n = 1 and —3 on large and small scales respectively, and with a 
turnover at some intermediate wavenumber that is determined by the single shape parameter 
F. Hereafter, we refer to this functional form of the power-spectrum as the "F model". We 
use it as a convenient parameterization, for comparability with other works, and for relaxing 
the COBE normalization. The free parameters that we vary in this case are the amplitude 
A and the shape parameter F. 



3.3. Error Models 

We make a special effort to estimate the distance errors. As mentioned in section |T|, 
this is done because the amplitude of the deduced PS depends on their sum in quadrature; 
if errors are overestimated, the PS is underestimated and vice versa. 

We first apply the likelihood analysis with the original distance errors, cXoj, as estimated 
a priori for each galaxy {i) in the SFI catalog with the procedure explained in section ^ 



Alternatively, we incorporate the errors in the likelihood analysis itself, by allowing a 
parametric model for the errors in addition to the parametric model of the PS. An error 
model is fully specified by the standard deviations cTj because we assume that the distance 
errors for the individual galaxies are uncorrelated and that the scatter is Gaussian. We try 
two alternative global modifications of the original errors as our error model: one is based 
on a free multiplicative factor, ctj = paoi, and the other is based on a free additive constant 



- 12 - 



in quadrature, cxj = (a^j ± q^y/"^. The latter is similar to the way we modeled the difference 
in scatter between the field and cluster samples (see § p.l| ). The errors are incorporated in 



the model that constitutes the correlation matrix, and the parameters p or q are adjusted 
simultaneously with the parameters of the PS until the likelihood is maximized. 

The apparent cost of adding the error parameter to the likelihood analysis is a larger 
formal error in the final results for the power spectrum and the cosmological parameters. 
However, since our original error estimates carries some uncertainty, this procedure, which 
provides an almost independent estimate of the errors, could add to the overall confidence 
in our results. 



4. TESTING THE METHOD 

4.1. Mock Catalogs 

It is essential to check the method with realistic mock catalogs, in view of the large 
errors in the data and the approximations made in the analysis. For this purpose, we use 
the A^-body simulation of Kolatt et al. (1996) which was designed to mimic the large-scale 
density distribution in our local universe. The simulation is based on initial conditions 
extracted from a reconstruction of the smoothed (5h^^Mpc Gaussian) real-space density 
field from the IRAS 1.2 Jy redshift survey, taken back into the linear regime. Small-scale 
perturbations were added by means of constrained random realizations, and the system was 
then evolved forward in time using a particle-mesh A^-body code until a present epoch defined 
by (Ts = 0.7. The "true" PS was calculated directly from the underlying mass distribution 
of the simulation, by Fourier transforming to k-space and calculating the power in bins of 
wavelength. 

"Galaxies" were identified in the simulation via a linear biasing scheme, and then divided 
into galaxy types, S's and E's, while obeying the morphology-density relation. Observational 
parameters were assigned to the S galaxies in the mock sample according to the prescription 
of Freudling et al. (1995), and perturbed at random according to the estimated observational 
errors. Subsequently, we selected ten random mock SFI samples using the exact selection 
criteria of the real SFI sample. Each of these mock catalogs was corrected for biases, and 
the errors were re-evaluated accordingly, in the same way as in the real data (see Eldar et 
al. 1998). 
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4.2. Testing with the T Model 

We first apply the likelihood analysis to the mock SFl catalogs using the F functional 
form (eq. [0]) as the prior model for the PS. We allow the amplitude A and the shape 
parameter F to vary, and include an additional free parameter in the error model. It is 
realized that the freedom provided by this family of models (just as by any other family of 
models) may not be enough for an adequate fit to the true PS. No additional constraint is 
applied on large scales, so this is a test of the ability of the velocity data alone to constrain 
the PS. 

Figure |1| (left panel) shows a contour plot of the resulting log-likelihood (ln£) in the 
parameter plane {A — T), as obtained from one of the realizations of the mock catalogs. 
The errors in this case were allowed to vary by the multiplicative factor p, and the plot 
shown corresponds to the best-fit error parameter. Here, and in all the figures that follow, 
the log-likelihood contours are relative to the maximum likelihood with contour spacing of 
A[ln£] = —1. The right panel of Figure ^ shows the corresponding best-fit power spectrum 
(solid line). The filled symbols mark the target of the reconstruction — the true PS of 
the simulation. The shaded area about the derived PS corresponds to the region of 90% 
confidence about the most-likely parameters in the likelihood plot, for fixed errors. The 
uncertainty becomes large at small fc's corresponding to scales larger than the sampled 
volume, because no additional data were used to constrain the PS on large scales. The 
figure demonstrates that for this random realization the likelihood analysis with the F model 
recovers the true PS well within the error-bars. A similar quality of recovery is obtained for 
all the random realizations of the mock catalogs, and also when the errors are varied in the 
alternative way. 

The maximum-likelihood errors are found to be within 5% of their "true" values. The 
latter were estimated by slightly modifying the known distance errors (as built into the mock 
catalogs) after correcting for Malmquist bias. The 5% error reflects the imperfect match 
between the assumed family of shapes for the PS and the true shape, and, perhaps, the 
uncertainty in the modification of the error estimate or the slight deviation of the modified 
errors from a Gaussian distribution. 



4.3. Testing with a Tilted ACDM Model 

We wish to check the success of the likelihood analysis also with the COBE-normalized 
CDM models. We choose as our test case the fiat {Q + = 1) ACDM family of models, 
with tensor fluctuations, a corresponding tilt in n, and a Hubble constant oih = 0.6. COBE 
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Fig. 1. — Left panel: Contour map of InC in the A — T plane for one random mock catalog. 
Contour spacing is A[ln£] = —1. A is in units of Aq = 2.0 x 10^( h^^Mpc)^. The errors were 
varied by a multiplicative factor. Right panel: The PS corresponding to the maximum-likelihood 
F-model parameters determined for the mock catalog. The filled symbols mark the true PS of the 
simulation. The shaded area around the PS is the 90% confidence region for the best-fit errors, 
obtained from the contour map. 

normalization is imposed as if the mock simulation is identical to the real universe. The 
likelihood analysis is thus performed by varying the parameters Q, n and the error-parameter 
p or q. This family of shapes for the PS is, again, not guaranteed to provide a perfect fit to 
the true PS. In particular, the parameter-dependent COBE normalization is not guaranteed 
to give the correct amplitude, since the simulation was not explicitly constrained to produce 
the level of large-scale CMB anisotropics detected in the real universe. 

Figure |^ (right panel) shows the best-fit power spectra of the 10 mock SFI catalogs, 
superimposed on the true PS. This test uses the q error parameter. The left panel shows InC 
contours in the Q—n plane for one representative mock catalog, with the maximum-likelihood 
points for all ten catalogs marked. We see that all the maximum- likelihood points fall along 
the ridge of high-likelihood in the one case plotted, and are therefore moderately consistent 
with one another. A way to translate the likelihood contours to errors in the values of the 
model parameters is by assuming that, with the errors fixed, — 21n£ has a distribution 
with two degrees of freedom. Then, the la confidence level around the maximum-likelihood 
point is at ln£ ~ —1.15 and the 90% confidence level is at ln£ ~ —2.3. The fact that indeed 
six of the ten cases fall within the la contour as determined above, and nine cases fall within 
the 90% confidence level, indicate that this crude error estimate is quite reasonable. The 
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Fig. 2. — Left panel: Contour plot of ln£ in the n — il. plane for one of the mock catalogs, for the 
tilted ACDM model with tensor fluctuations and h = 0.6, and the errors varied in quadrature. The 
best-fit values for all catalogs are marked by '+'. Right panel: Best-fit PS of the 10 mock catalogs 
(thick lines representing same curves derived for two different catalogs) . The shaded area represents 
the 90% confidence region for the catalog whose contour plot is shown. The filled symbols mark 
the true PS of the simulation. 



90% confidence region for this specific catalog is again drawn as a shaded area in the PS 
plot; one can see that this region indeed resembles the actual scatter of the ten cases. 

The maximum-likelihood power-spectra fit reasonably well the true PS, with a fairly 
small spread on small and intermediate scales. For large scales (small fc's) the scatter is 
somewhat larger, but not as large as for the P model which was completely free at large 
scales. Again, the success of recovery is similar when the alternative error parameter is used, 
and the errors are similar to those obtained in the case of the P model. It is encouraging 
to note that the recovery of the PS is fairly robust on the relevant scales {k ~ 0.1 /iMpc~^) 
among the realizations, and independent of the prior model assumed for the PS, or the 
assumed error model. 



5. RESULTS 



5.1. Maximum-Likelihood Errors 



Before estimating the PS from the actual SFI data, we investigate the reliability of our 
observational error estimate by allowing certain freedom in the errors. As a test case we use 
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Fig. 3. — Left panel: Contour plot of InC in the 0, — n plane for the SFI sample, for the 
tilted ACDM model with tensor fluctuations and h = 0.6, using the original error estimates. The 
best-fit point is marked with a '+'. The maximum likelihood locations when varying the errors in 
quadrature ('g') or by a multiplicative factor ('p') are also marked. Right panel: The most-likely PS 
for this model, for these three variants of the errors. The dotted region around the PS represents the 
90% confidence limit for the case of the original errors (thick line), obtained from the high-likelihood 
ridge shown in the contour map. 

as a prior for the PS the COBE-normahzed ACDM family of models, with tensor fluctuations 
and a corresponding tilt in n, and with the Hubble constant fixed at h = 0.6. We perform 
the likelihood analysis on the real SFI data varying Q and n, with the errors treated in three 
different ways; first with the errors fixed at their original values, and then by varying them 
according to the two error models discussed above. 

Figure ^ summarizes the results obtained in these cases. The left panel shows the 
In— likelihood contours in the Q — n plane for the case of fixed errors, with the best-fit points 
for the three cases marked. The corresponding power spectra are presented in the right 
panel. In the p case, the preferred errors are 5% larger than the original ones, while in the 
q case the preferred errors are smaller by 0.03 in quadrature (typically a decrease of ~ 2%). 

The different trends in the likely errors refiect our uncertainty of the exact form of the 
error model. We note that while these changes are in opposite directions, they are of small 
magnitude, within the uncertainty expected based on the mock catalogs. The corresponding 
changes in the best-fit parameters are along the ridge of high likelihood in the Q — n plane, 
within the la confidence level, i.e., it is hardly significant. In all three cases, x^/^dof ~ 1 for 
the best-fit PS (1.02 for the original, fixed errors, 0.99 for the p error model and 1.02 for the 
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q error model), implying that all are reasonable fits to the data. Similar results concerning 
the errors are obtained when the other PS models are used as priors. The error estimate is 
robust to variations in the original errors about which the error model is perturbed. This 
likelihood analysis of the errors thus provides a very encouraging indication that the original 
error estimates in SFI are accurate to better than 5%. Note that "original" here refers to the 
refined SFI errors after the correction for biases (§ p.2|) . The fact that the likelihood analysis 
and the semi-analytic correction converge to the same error estimate is encouraging. Based 
on this finding, we perform the rest of the analysis in this paper using fixed errors at their 
original values. 



5.2. COBE Normalized CDM Models 

We now use the generalized CDM families of cosmological models of the form described 



in equations @ and ([T0|). Our models include open CDM (OCDM ), fiat models with a 
cosmological constant, and tilted models with or without tensor fluctuations, allowing for 
variations in the cosmological parameters Q, h, and n. For each specific choice of model and 
parameters the amplitude is fixed according to the 4-year COBE normalization. 



5.2.1. Scale- Invariant Models 

Figure |^ shows the resulting likelihood contours for the scale- invariant case, n = 1, for 
the OCDM model and the ACDM model. The contours are plotted in the Q ~ h plane. 
The best-fit parameters in each case are marked, but it is clear from the elongated contours 
that the two parameters are not determined separately. The high-likelihood ridges rather 
constrain a degenerate combination of these parameters, which can be roughly fitted by the 
following functions: 

^^/iggO-s = 0.68 ±0.06, OCDM; (13) 
nheo^-^ = 0.59 ±0.07, ACDM. (14) 

The error-bars (here, and throughout the paper) arise from the joint 90% confidence region of 
the parameters. The constraints on asf{^), obtained by integrating over the corresponding 
power spectra, are O.SS^'^q;?]' and 0.8llo;j7 for OCDM and ACDM respectively. The error-bars 
quoted are the marginalized 1-dimensional 90% confidence limits. For an assumed value 
of h, e.g., h = 0.6, the maximum-likelihood values are Q = 0.67 ± 0.05 for OCDM and 
n = 0.58 ± 0.06 for ACDM (as marked on the plots). The xV^dof for the best-fit PS 
are 1.01 and 1.04 respectively, with similar values along the high-likelihood ridge. With 
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Fig. 4. — Contour plot of ln£ in the ^l — h plane for scale-invariant OCDM model (left panel) and 
ACDM model (right panel). The most-likely value of 17 and its 90% error bar are marked for a 
fixed value of h = 0.6. 



-^dof = 1213, one expects for a good fit 1.00 ±0.04, so our CDM models indeed provide good 
fits to the data. 

5.2.2. Tilted Models 

Figure |^ presents the results obtained when allowing for a tilt in the PS on large scale 
relative to n = 1, both for the OCDM and ACDM families of models. The first cases 
considered are with scalar fiuctuations only, T/S = 0. We fix the Hubble constant here at 
h = 0.6 while varying Q and n. Again, the elongated ridge of high-likelihood determines a 
certain degenerate combination of the parameters, which can be approximated by: 

fin^-^(/i6o°-^) = 0.68 ±0.07, OCDM; (15) 
nn'^-^ (heo^-^) = 0.58 ±0.08, ACDM. (16) 

The h dependence is determined for the n = 1 case. The corresponding constraints are 
asf{n) = 0.83lHo for the tilted OCDM case and asf{n) = 0.82+1^^ for tilted ACDM . The 
X^/A^dof values are 1.02 in both cases, again a good fit. 



The case of the tilted ACDM family of models, with h = 0.6 and with a tensor 
component of T/S = 7(1 — n), was partly discussed already as our default case in § pTl . 
The likelihood map in Figure ^ reveals the familiar situation of a high-likelihood ridge that 
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11 n 
Fig. 5. — Contour plot of InC in the — n plane, for the tilted OCDM model (left panel) and 
tilted ACDM model (right panel). In both cases h = 0.6 and no tensor component is included. 



constrains a degenerate combination of the cosmological parameters, now approximated by 

Qn^-^iheo^-^) = 0.58 ± 0.08 , ACDM + tensor . (17) 

The h dependence is determined for n = 1. The corresponding value of asfi^^) is 0.8lto;o8. 
The uncertainty associated with the PS, shown as the shaded area in the right panel of 
Fig. ^, is similar to the uncertainty in the other COBE-normalized CDM variants. 



5.3. r Model 

Finally, we use the F model as a prior for the PS, varying the amplitude A and shape 
parameter F with no additional constraints imposed at large scales. Figure ^ shows the 
contours of ln£ in the A — T plane, and the corresponding best-fit PS. The maximum 
likelihood values are F = 0.375 ± 0.14 and A = 5.0 x 10^( h~^Mpc)^. The P^i' degree 
of freedom for the maximum likelihood parameters is x^/^dof = 1-03, indicating that 
this is a good fit to the data. The constraint obtained by integrating over the power 
spectra is (j8/(f2) = 0.80l'^o;o8. The scatter at small /c's is larger than in the case of the 
COBE-normalized models, due to the amplitude freedom, as seen already in the mock 
catalogs. 
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Fig. 6. — Left panel: Contour plot of InC for the T model. The best-fit point is marked with a 
'+'. Right panel: The best-fit PS, with the shaded area marking the uncertainty based on the 90% 
confidence region of the likelihood contours. 



6. ROBUSTNESS OF RESULTS 



The error estimates in the parameters given in the previous section are formal 90% 
confidence levels. In this section, we test the robustness of these results to various variations 
in the data and models used. 



6.1. Robustness to Models 

Figure ^ shows the power spectra corresponding to the maximum-likelihood parameters 
for all the models presented so far in this paper, including the COBE-normalized CDM 
variants and the F model. The 90% confidence region for the tilted ACDM model with 
tensor fluctuation and h = 0.6 is drawn as well, as a reference for the uncertainty associated 
with each model based on the likelihood contours. The similarity of all the curves is striking; 
they agree well within the formal uncertainties of each other. The agreement is excellent 
for k > 0.1, on the scale where the data constrain the models effectively. The difference 
between the curves shows a slightly larger scatter on larger scales, not properly sampled by 
the present data. The similarity of the results using as priors the COBE-normalized CDM 
models and the amplitude-free F model indicates that the peculiar velocity data themselves 
contain meaningful information to constrain the PS. 
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Fig. 7. — The maximum-likelihood power spectra based on the F model (solid line) and the various 
COBE-normalized CDM models. The shaded area is once again the 90% confidence region for the 
tilted ACDM model with tensor fluctuations and h = 0.6. 



Table 1 summarizes the features of the most likely power spectra based on the various 
prior models. The approximate constraint on the combination of cosmological parameters 
as obtained from the high-likelihood ridge is given for each case. The best-fit values of 
the individual cosmological parameters are also listed, but recall that they carry large 
uncertainties. The exact location of the maximum-likelihood point in the high-likelihood 
ridge is hardly significant. Parameters that were held fixed in the likelihood analysis are 
marked in brackets. Several characteristics of the best-fit power spectra are listed: the value 
of asf{^), the amplitude of f{fiyP{k) at k = 0.1h~^Mpc and the location of the PS peak, 
kpeak- The error-bars quoted in the header represent the typical 90% confidence uncertainty 
in these quantities within each family of models. 

The typical results for the PS are P{k = 0.1 /iMpc-^)^^ ^ = (4.4 ±1.5) x 10^ {h-^Mpcf 
and agfl^'^ = 0.82 ± 0.10. The variations from model to model are much smaller than 
the formal errors for each model, increasing the above errors to 1.7 and 0.12 respectively. 
These results are thus almost independent of the model, at least for the family of models 
considered here. The actual likelihood values of all the best-fit models are very similar, and 
all have comparable x^/^dof — 1 values. The variation of the high-likelihood ridge between 
the ACDM and OCDM families of models is more noticeable. The general constraint on 
the combination of cosmological parameters can be roughly approximated by Qn" h^Q^ = 
0.62 ±0.15, where the error includes the formal uncertainties of the three parameters and the 
variations between models. For ACDM , fi = 1.3 and u = 2.0, 3.9 without and with tensor 
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Table 1: Maximum-Likelihood Results for the various models 

COBE-normalized CDM models 











^peak 














High Likelihood 




{h-^ Mpc^){hMpc-^ 


) 










CDM Model 


Ridge 


(±0.10) 


(±1500) 


(±0.01) 


n 


n 


h 




-ln£ 


Open, n=l 


VLhQQ^-^ = 0.68 ± 0.06 


0.83 


4400 


0.038 


0.55 


(1) 


0.78 


1.01 


8579.4 


A, n=l 


nheo^-'^ = 0.59 ± 0.07 


0.81 


4600 


0.031 


0.95 


(1) 


0.42 


1.04 


8580.1 


Tilted-Open 


Qn^-^ = 0.68 ±0.07 


0.83 


4600 


0.035 


0.72 


0.96 


(0.6) 


1.02 


8579.5 


Tilted-A 


Qn^-o = 0.58 ±0.08 


0.82 


4200 


0.037 


1.00 


0.76 


(0.6) 


1.02 


8579.6 


Tilted-A 


Qn^-^ = 0.58 ±0.08 


0.81 


4300 


0.037 


0.79 


0.92 


(0.6) 


1.02 


8579.5 


+tensor 




















r model 
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fluctuation respectively. For OCDM , without tensor fluctuations, the powers are fi = 0.9 
and u = 1.4. 

The similarity of the power spectra obtained using the COBE-normalized CDM models 
and the COBE-free F model (see Table 1) indicates that the PS is predominantly determined 
by the velocity data. Therefore, we have so far ignored the error associated with the COBE 
normalization. As a test for the sensitivity to this error, we have repeated the analysis using 
the tilted ACDM model (with tensor fluctuations), but now normalized alternatively ~ 18% 
higher or lower then the mean COBE values (in accordance with the relative ±la uncertainty 
associated with Qrms-ps\n=i, Bennett et al. 1996). This results in a slight shift of the 
high-likelihood ridge, corresponding to a ~ 6% change in the constraint on the combination 
of parameters (eq. |17]; a smaller value is obtained for the higher normalization and vice 



versa), which is within our formal Icr error-bars. However, the combined effect of the different 
amplitude and corresponding cosmological parameters on the PS is essentially negligible, 
with CTg^^'^ varying by only 0.01. 



6.2. Zero-point Uncertainty 



A fundamental freedom in the measured peculiar velocities is in the global zero-point 
of the TF relation, which fixes the distances at absolute values (in kms^^). Changing the 
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Fig. 8. — Left panel: Contour plots in the Q — n plane for the tilted ACDM h = 0.6 model 
with a tensor component, for the original calibration. The maximum-likelihood point is marked 
by '+'. The maximum-likelihood values when varying the global zero-point by a ±5% Hubble flow 
are marked by '-I-5V-5' accordingly. Right panel: The most likely PS for the original zero-point 
calibration (solid line) and when varying the zero-point by e = -1-0.05 (long dashed line) and —0.05 
(short dashed line). The corresponding agCl^'^ values are marked on the plot. The shaded region 
is the formal likelihood 90% confidence region for the original case. 

zero-point, that is multiplying the distances r by a factor (1 — e), is equivalent to adding a 
monopole Hubble-Iike flow er to the peculiar velocities. The zero-point calibration of the TF 
relation used for the SFI sample was obtained from the SCI catalog of ~ 500 galaxies within 
24 clusters, using the "Basket of Clusters" approach (Giovanelli et al. 1997a, 1997b). The 
uncertainty in the zero point was estimated to be about 0.05 magnitudes, which corresponds 
to an uncertainty in the velocity field of 2.5% of the distance. 

To estimate the effects of such uncertainties we have run the likelihood analysis with 
our tilted ACDM test-case, conservatively using zero-point changes of twice the estimated 
uncertainty, e = ±0.05. Figure |^ illustrates the effect on the results for these cases. 
The changes of zero-point appear to shift the location of the maximum-likelihood values 
essentially along the high-likelihood ridge. The high ridge is not altered by much when the 
zero point varies in this range. The right panel shows the resulting best-fit PS for the three 
different zero points, and lists the corresponding values of a^^*^'^ ■ While the variations in 
the zero-point systematically affect the PS, the changes are not large; they fall within the 
range of the formal likelihood errors, and are of the same order of the uncertainty associated 
with the random distance errors (compare to Fig. |^). It is encouraging that the amplitude 
of the PS on intermediate scales {k ~ 0.1 /iMpc"^) is robust vis-a-vis changes in in the zero 
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point. Similar results were obtained when using the other families of PS prior models. 

Similar to the uncertainty in the zero point of the distance indicator, there is also an 
uncertainy associated with the slope of the TF relation. This could lead to correlated errors 
in the inferred distances and peculiar velocities, due to the fact that the average linewidth 
of SFI galaxies slightly depends on distance (Wegner et al. 1999). However, the impact of 
this uncertainty on our results is even slightly smaller than that of the uncertainty in the 
zero point, perhaps because the SFI sample has been selected intentionally to minimize the 
distance dependence of the line widths. 

6.3. Nonlinear Effects 

A basic assumption in our analysis has been that linear gravitational instability theory 
is adequate for the purpose of recovering the PS from observed velocities on the scales of 
interest here. This is based on the fact that in the mildly-nonlinear regime the velocity field 
is approximated by linear theory better than the density field (basically because the velocity 
is a spatial integral of the density and is affected by fluctuations on larger scales). Indeed, 
the success of the recovery of the PS from the mock catalogs leads us to believe that this 
assumption is justified. However, one cannot rule out the possibility that some nonlinear 
effects are artificially reduced to some degree in the particle-mesh A^-body simulation, and 
it is possible that the smooth shape of the linear PS as predicted for the CDM family of 
models may fail to properly match the nonlinear features that may be present on small scales 
in the real data. Therefore, we discuss in this section possible nonlinear effects, which could 
manifest themselves in different forms. For example, as coherent motions associated with 
the non-linear evolution of the PS (§ |6.3.1| ), or as incoherent random motions, perhaps due 
to shell-crossing, which may be modeled as an additional velocity component of dispersion 

6.3.1. Nonlinear Power Spectra 

A way to include more properly nonlinear effects in our analysis is by developing 
an approximation for the nonlinear evolution of the PS and then incorporating it in the 
likelihood analysis. Such approximations exist for the density power spectrum, Ps {e.g., 
Peacock & Dodds 1994; Jain, Mo & White 1995; Peacock & Dodds 1996, hereafter PD), 
but we need a similar approximation for the evolution of the velocity power spectrum, P^,, 
which is the quantity we actually confront with the data. A development and application of 
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such an approximation is beyond the scope of the present paper and will be presented later 
(Zehavi et al. 1999). Here, we summarize some relevant issues and illustrate the magnitude 
of such effects. 

Figure |^ shows the velocity PS computed in several different ways from an adaptive 
P^M cosmological A^-body simulation with a resolution higher by an order of magnitude 
than the simulation used in the present paper for the mock catalogs (but inside a smaller 
box of size 85h~^Mpc; GIF simulation, Colberg et al. 1999). The initial model used for is 
the so-called rCDM model with Vt = 1.0, h = 0.5 and a modified shape parameter P = 0.21. 
The figure clearly demonstrates that the velocity PS is reproduced by linear theory much 
better than the density PS. The Py that is computed directly from the evolved velocity field 
of the simulation (solid dots) lies slightly below the Py obtained from the assumed Ps using 
linear theory {Py oc k~'^Ps, solid line). On the other hand, the nonlinear correction to Ps 
{e.g., PD) is larger than that of Py and in the opposite direction (upwards, as can be seen 
by the open dots and dashed line in Fig. 

One might have naively expected that the likelihood analysis using a pure linear 
treatment would be inferior to incorporating a non-linear correction for Ps followed by a 
linear translation to Py. However, as illustrated in Figure ^, this is not the case. The latter 
procedure overestimates the nonlinear effects on the velocity PS and increases the bias in the 
results. A similar bias is reproduced when using the mock catalogs from the low-resolution 
simulation of Kolatt et al. (1996), which exhibits a similar behavior as in Fig. ^. This could 
be remedied, in principle, by incorporating the evolved Ps and then counter-balancing it 
with a proper approximation for the nonlinear velocity-density relation, but this would be 
risky as we would be applying two large corrections in opposite directions to mimic a small 
net effect. Until we develop a direct nonlinear correction for Py, we adopt the fully linear 
procedure as our best approximation. This is justified by its success in the mock catalogs 
and by the expectation for only small nonlinear effects in Py. 

6.3.2. Random Motions 

We have made an ad-hoc attempt to model non-linearities by introducing an 
uncorrelated velocity component of constant dispersion ay, that adds a free term at zero-lag 
to the correlation function derived from the linear PS model. This may be a crude way 
to represent small-scale random motions that are associated with multi-streaming. An 
alternative interpretation of this additional parameter may be as an unrecognized uncertainty 
in the distance estimate which does not depend on distance and is therefore not included in 
our usual error model. In either case, this provides a test for the robustness of our results 
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Fig. 9. — The velocity PS as computed from tlie liigh-resolution simulation of Colberg et al. (1999) 
(solid dots), compared with the theoretical linear PS (solid line), the corrected PS using the PD 
formalism (dashed line) and a computation via the density PS of the simulation (open dots). The 
latter three have been transformed to velocity PS using the linear velocity-density relation. 



to an additional degree of freedom. 



Figure |10| demonstrates the effect of including a free in the likelihood analysis, again 
for our tilted ACDM test-case. When allowing for this extra freedom, the preferred value 
turns out to be cr^ = 200 ± 120 kms~^, and is associated with a PS that is slightly lower 
for k > 0.1 and somewhat higher at small k. The value of ag^^'^ is reduced by 14%. The 
deviations, in general, are comparable to the formal likelihood 90% uncertainty marked by 
the shaded area. The likelihood contours are somewhat sparser in this case, because of the 
additional scatter that reduces the sensitivity to variations in the parameters. The ridge of 
high likelihood is slightly shifted toward smaller values of the cosmological parameters and 
it can now be roughly described by 

nn^-^iheo^-^) = 0.49 ±0.09, (18) 

a ~ 15% decrease compared to equation (p!?]), which is of the order of the random error. 



Such a preference for a non-zero cr„ associated with a change in the PS is not recovered in 
the mock SFI catalogs, for which a similar likelihood analysis turns out to prefer a negligible 
(T„ and a negligible effect on the PS. A^-body simulations of higher resolution may clarify 
this situation. 

We note that the inclusion of a free ay in the fit to the real SFI data leads to results 
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Fig. 10. — Left panel: Likelihood contour plot in the 0— n plane for the tilted ACDM h = 0.6 model 
with tensor fluctuations, with an additional scatter of = 200 kms~^. The maximum-likelihood 
point is marked by '*', and the corresponding point for cjt, = is marked by '+'. Right panel: The 
most likely PS when including o"^ = 200 in the fit (dashed line) and for the original cr^ = case 
(solid line) together with its 90% confidence region shaded. 

similar to those obtained when including a free multiplicative parameter in the error model 
(§ ^.1|) . The interpretation of a nonzero is thus not unique: It may refer to nonlinear 
effects that exist in the real data but not in the current simulation, or it may indicate that 
the actual errors are slightly larger than the original estimates. Since there is no clear benefit 
from adding this extra parameter and the theoretical justification as a model for nonlinear 
effects is weak, its inclusion in our main-stream analysis does not seem to be justified. Still, 
in our total error-balance, we consider a systematic error of 15% due to non-linear effects. 



6.4. Comparison to the PS from Mark III 

A similar likelihood analysis (though with errors fixed a priori) has been recently applied 
by Zaroubi et al. (1997) to the Mark III catalog of peculiar velocities. (Willick et al. 
1995; 1996; 1997a). It is interesting to investigate whether the recovered power spectra 
are consistent with each other, given the respective uncertainties. This is intriguing because 
there are certain differences in the velocity fields as reconstructed from the two samples, 
especially in the bulk flows both in the very local neighborhood and of outer shells {e.g., da 
Costa et al. 1996; 1998; Dekel 1998; Dekel et al. 1998; Giovanelli et al. 1998a, 1998b). 
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Fig. 11. — Comparison of SFI PS results to Mark III for the tilted ACDM test-case. Top left panel 
shows ln£ contours for the SFI data, top right panel shows the contours for the Mark III data. The 
best-fit parameters for SFI and Mark III are marked, on both, by 'S' and 'M' respectively. The 
lower panel shows the maximum likelihood PS corresponding to SFI (solid) and Mark III (dashed). 
The 3 solid dots mark the PS calculated from Mark III by Kolatt and Dekel (1997), together with 
their estimated la error-bar. The shaded region is the SFI likelihood 90% confidence region. 
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Figure 11 presents a comparison of the likelihood analysis for the two samples using 
our representative tilted ACDM family of models. Shown once again are the likelihood 
contours for the SFI data, together with the corresponding plot for the Mark III data. The 
high-likelihood ridge is similar for both samples. While the Mark III result slightly favors 
higher values of fl and lower values of n, the differences are along the ridge of maximum 
likelihood and are therefore hardly significant. The contours are slightly more concentrated 
for the Mark III catalog because it consists of more galaxies. 



The best-fit PS for the two catalogs are shown in the bottom panel of Figure [TJ, on 
top of the shaded area which marks the 90% confidence region for SFI. The resultant power 
spectra are consistent within the errors, and they agree particularly well on intermediate 
scales, where the data provides the most meaningful constraints. The corresponding best 
values for agQ'^'^ are 0.81 and 0.85 for SFI and Mark III respectively. Similar results are 
obtained when comparing likelihood analysis of the two catalogs using the other PS models. 
It is worth noting here that the systematics discussed in the previous sections, with regard to 
the SFI analysis, are found to affect the Mark III likelihood analysis in a similar way. The PS 
computed by Kolatt & Dekel (1997) from the Mark III smoothed density field recovered by 
POTENT is also displayed on the figure (as three symbols with error bars). The agreement 
of the SFI result with this independent calculation of the Mark III PS is good. 

A recent comparison of Mark III with IRAS 1.2 Jy (Willick & Strauss 1998) suggests 
an alternative zero-point calibration for one of the Mark III datasets. We have applied our 
likelihood analysis to the Mark III data revised accordingly, and found negligible changes in 
the resulting power spectrum and cosmological parameters, smaller than the uncertainties 



due to global zero-point discussed in § 3.2 



The close agreement between the mass power spectra derived from the two datasets 
indicates that the results presented here are quite robust and are unlikely to arise form 
specific peculiarities of either of the two samples. This does not preclude possible differences 
that are not picked up by the specific statistic used - in our case, the mass PS. In particular, 
the difference in the two bulk flows, which is known to exist, is not reflected in the power 
spectra. This is because the wavenumbers corresponding to the bulk velocity are smaller 
than the k range that dominates the fit in our current analysis. 



7. CONCLUSION 



We used a linear maximum-likelihood method to measure the mass-density power 
spectrum from the SFI catalog of peculiar velocities, and to determine the cosmological 
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parameters for families of physical CDM models with or without COBE normahzation. We 
have corrected for biases introduced by the non-trivial selection procedure of the SFI catalog 
using a new semi-analytic procedure. We have verified that the results are quite insensitive 
to the detailed way by which we implement this bias correction. This approach allows also 
to refine the distance errors estimates. Our new version of likelihood analysis enabled us to 
independently verify the error estimates of the SFI catalog to within an uncertainty as small 
as 5% of the error, which we regard as very encouraging. Since the errors affect the PS in 
a systematic way, this independent confirmation adds significantly to our confidence in the 
results. 

The general result for all the models examined here is that the power spectrum at 
k = 0.1/iMpc-i is P{k)n^-^ = (4.4 ± 1.7) x 10^ {h-^Mpc)\ and that ag^^'^ = 0.82 ± 0.12. 
These results are obtained by the peculiar-velocity data independent of the specific shape 
assumed for the PS, and are consistent with the result of the F model independent of the 
COBE normalization. The random errors quoted are 90% confidence level and they include 
small variations due to the choice of model for the PS within the families of models tried 
here. 

For the general family of COBE-normalized CDM models, we find a high-likelihood 

ridge in the Q ~ n — h parameter space, which can be crudely approximated by fln'^ h^o^ ~ 
0.62 ± 0.15, where for ACDM n = 1.3 and u = 2.0, 3.9, without and with tensor fluctuation 
respectively. For OCDM , without tensor fluctuations, the powers are fi = 0.9 and u = 1.4. 
Again, the error quoted is the formal 90% uncertainty including the model variations. Thus, 
for h = 0.6, the maximum-likelihood value of fl ranges between 0.6 and unity while n varies 
between 1 and 0.8 respectively. Without a tilt, values of fl as low as 0.5 are allowed within 
the 90% confidence limit. 

Our tests using mock catalogs based on an iV-body simulation that mimics our 
cosmological neighborhood indicate that the systematic errors in our results are relatively 
small. In particular, the nonlinear effects in the mock catalogs are found to be negligible. 
This is indeed expected because the quantity we actually measure is the velocity power 
spectrum in the mildly-nonlinear regime, which we have demonstrated to be reasonably 
approximated by linear theory. An ad-hoc test for nonlinear (multi-streaming) effects in 
the data themselves indicated that they may work to reduce the values of the cosmological 
parameters given above, but that this effect is not larger than ~ 15%. In order to reflne our 
estimates of the systematic effects even further, we intend to repeat the current analysis using 
a proper nonlinear scheme, and to repeat the tests of the method using simulations of higher 
resolution which are in preparation. We thus estimate the total systematic uncertainty to 
be of order ~ 15%, namely comparable in size to the random errors. Therefore, to be on 
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the safe side when comparing our results to other results, we recommend as a rule of thumb 
multiplying the quoted errors by a factor of ~ 2. 

As yet another word of caution, it is worth recalling that our analysis is heavily weighted 
by the galaxies at relatively small distances, because the data is weighted by the inverse 
squared of the distance errors. This means that the result is sensitive to the data and 
error estimate of the inner galaxies. It is possible in principle that a source of distance 
error which operates preferentially at small distances has somehow escaped our attention 
and is not properly modeled by our error model. To test the effect of such a possibility, 
we have repeated the analysis after pruning all galaxies with distances smaller than a given 
distance. When pruning inside 15 h^^Mpc (3% of the data) we obtain for the most likely 
value CTgQ^-^ = 0.85 instead of the original result of o-gJl^-^ = 0.81 when using all the data (still 
with our standard pruning based on linewidth). When pruning inside 25h~^Mpc (17% of 
the galaxies), we obtain instead a^O,'^'^ — 0.71. It is encouraging to find that these variations 
are within the 90% hkelihood contours of the different cases in the Q — n plane, but this is 
yet another potential source of uncertainty to bear in mind. 

A systematic trend does seem to show up when we eliminate as much as the whole 
inner half of the data (inside a distance of 46h~^Mpc, or with hnewidth smaller than 2.48); 
the outer data, when analyzed by themselves, indicate a significantly lower PS then the 
inner data. This effect is not reproduced in the mock catalogs and is therefore not likely 
to represent a general fault in the method. Possible explanations for this effect are larger 
uncertainty in our estimate of random and systematic errors at large distances, differences 
between the assumed and the true TF relation, or a true PS with a different shape than 
our models. It may also be due to a real difference between the density fields in the two 
halves (that is somehow not properly reproduced in the simulation) or to a systematic 
dependence of velocity bias on galaxy properties. We carried out a number of tests in which 
we added to the likelihood analysis ad-hoc free parameters which allow more fiexibility in 
the distance dependence. These include variations in the TF parameters and in the errors 
as a function of linewidth. Our tests indeed led to some improvement in the agreement 
between the two subsamples, but, being only preliminary, they did not yield so far a firm 
conclusion as for the dominant source of the effect and the optimal way to deal with it. 
Since the variations introduced preferentially affect the peculiar velocities of galaxies at 
large distances, which typically have large errors and therefore contribute only little to the 
likelihood procedure, they do not affect significantly the resultant power spectrum from the 
full sample. We therefore conclude that our current results are robust, and defer a more 
thorough investigation of this trend to a future analysis. 

The recovered mass power spectrum, and the constraints on the cosmological parameters 
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obtained here, are consistent with the results of a similar analysis apphed to the Mark III 
catalog of peculiar velocity. This is despite the fact that these two catalogs seem to differ in 
some of their other properties, such as the large-scale bulk velocity. Indeed, the bulk velocity 
is not expected to contribute to the density on smaller scales. There is also an apparent 
disagreement between the results obtained from peculiar velocities of clusters (Borgani et al. 
1997) and our result for the SFI field galaxies. 

As mentioned in the Introduction, our dynamical result of erg = a^Vf*'^ ~ 0.8 ± 0.2 
may be crudely compared to estimates of the j3 parameter obtained when comparing the 
same SFI data to a rcdshift survey of galaxies, da Costa et al. (1998) find (3 = 0.6 ± 0.1 
when comparing the SFI peculiar velocities to the velocities predicted by the IRAS 1.2 Jy 
redshift survey, assuming linear biasing. A similar value was obtained from Mark III when 
the comparison was done via velocities (Davis et al. 1996). With asg ^ 0.7 for IRAS galaxies, 
the predicted /3 from our current constraint on erg, via /3 = ^s/'^sgi is significantly closer to 
unity than to 0.6 (compare also to Kolatt & Dekel 1997, Fig. 6). The residuals between 
the measured peculiar velocities and the IRAS predictions, for the best-fit (3 value, were 
found in this comparison (da Costa et al. 1998) to be significantly higher than the errors 
as originally estimated for the SFI data (based on the scatter observed in the SCI cluster 
sample) combined with the errors estimated for the IRAS data. 

One possibility is that the the IRAS model fails to predict some of the peculiar 
velocities that exist in the SFI data, e.g., because the distribution of galaxies is not properly 
approximated by a simple, linear, scale-independent and deterministic biasing relation {e.g., 
Dekel & Lahav 1998). In that case, the interpretation of the value of j3 determined from 
fitting the IRAS predictions to the velocity data is not clear. 

On the other hand, it is also possible that the errors in SFI arc indeed larger than 
originally estimated. Such larger errors would accordingly reduce the PS amplitude, in 
particular the value of ag, as estimated in the current paper. However, such an effect should 
have been detected by investigating the global biasing properties of the sample. It would 
also be hard to understand why our likelihood analysis does prefer errors very similar to 
the original estimates. Although we are fairly convinced that nonlinear effects in the current 
analysis are confined to the level of < 15%, it will be worth making an extra effort to improve 
the accuracy in a future paper. 
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